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We propose angle-resolved two-photon photoemission spectroscopy (AR-2PPES) as a technique to 
detect the location of the bottom of the upper Hubbard band (UHB) in two-dimensional insulating 
cuprates. The AR-2PPES spectra are numerically calculated for small Hubbard clusters. When 
the pump photon excites an electron from the lower Hubbard band, the bottom of the UHB is less 
clear, but when an electron in the nonbonding oxygen band is excited, the bottom of the UHB can 
be identified clearly, accompanied with additional spectra originated from the spin-wave excitation 
at half filling. 

PACS numbers: 71.10.Fd, 79.60.-i, 74.25.Gz 



The charge gap in Mott insulators is a consequence 
of strong electron correlation. The nature of excitations 
across the gap is controlled by both the lower Hubbard 
band (LHB) and the upper Hubbard band (UHB). There- 
fore, the clarification of the momentum dependence of 
the two bands will be very crucial for understanding the 
nature of the Mott-gap excitation. 

Parent compounds of high-T c superconductors such as 
La2Cu04 and Ca2Cu02Cl2 are a good example of the 
Mott insulator in two dimensions (2D). It was theoreti- 
cally proposed that the top of the LHB [more precisely 
the Zhang-Rice singlet band (ZRB) (Ref. 1)] is located 
at k = (±tt/2, ±tt/2), while the bottom of the UHB is at 
[ir, 0) and (0, 7r)£ The location of the top of the LHB has 
clearly been observed by angle-resolved photoemission 
experiments^ On the other hand, the bottom of UHB 
has not been directly observed yet, but only indirectly 
confirmed by resonant inelastic x-ray scattering that re- 
veals momentum-dependent Mott gap excitations 4> 

In this Brief Report, we propose angle-resolved two- 
photon photoemission spectroscopy (AR-2PPES) as a 
new technique to detect the location of the bottom of 
the UHB in the 2D insulating cuprates. The two-photon 
photoemission spectroscopy has been widely used for the 
studies of the electronic excitations at metal surfaces and 
their decay in the time domain^^ On the other hand, 
for strongly correlated systems such as high-T c cuprates 
there are a few studies only in the metallic and supercon- 
ducting regions of Bi 2 Sr2CaCu208+5r i ^ but no report on 
insulating materials. Therefore, the application of AR- 
2PPES to the insulating cuprates will open a new way of 
spectroscopic study for strongly correlated electron sys- 
tems. 

The AR-2PPES spectrum is numerically calculated for 
small clusters described by the Hubbard model with real- 
istic parameters. We consider two types of pump-photon 
excitation: One corresponds to the case where the pump 
photon induces an excitation from LHB to UHB, while 
the other is an excitation from the nonbonding (NB) oxy- 
gen band in which the oxygen 2p a orbitals form an asym- 
metric state with respect to the Cu3d x a_y2 orbital. In the 
former, we find two kinds of AR-2PPES spectra that re- 



flect information on either the UHB or LHB. However, 
the location of the bottom of UHB is not clearly identi- 
fied because of diffusive features in the spectra. In the 
latter case, we clearly find the location of the bottom of 
UHB. In addition to this, new spectral weights that are 
not seen in inverse photoemission spectra emerge in AR- 
2PPES. The origin of the additional spectra is attributed 
to the spin- wave excitation expected at half filling. These 
theoretical results will be useful for the analysis of AR- 
2PPES experimental data in the near future. 

We employ a single-band Hubbard model with long- 
range hoppings to describe the electronic states in the 
2D insulating cuprates. The Hamiltonian reads 
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where c\ a is the creation operator of an electron with spin 

cr at site i, and n\^ a — c\ _Ci )Cr . The summations run 
over neighboring pairs up to the third-nearest neighbors. 
The hopping parameter Uj then contains three terms: t, 
t', and t" for the first, second, and third neighbors. U 
is the on-site Coulomb interaction, and fi is the chemical 
potential. We take the parameter values of the t-t'-t"-U 
model to be t'/t = -0.34, t"/t = 0.23, and U/t = 10, 
which are realistic ones for the 2D cuprates^ 

In AR-2PPES examined in the present work, we con- 
sider the situation that the pump photon induces a dipolc 
transition from occupied states to UHB and the probe 
photon kicks an electron in UHB out of sample. Gen- 
erally there are two kinds of contributions for the elec- 
tron energy distribution of AR-2PPES>2iiS One is due to 
simultaneous excitations that is given, at zero tempera- 
ture, by 
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FIG. 1: Schematic picture of two-photon photoemission pro- 
cesses in insulating cuprates. Two processes are shown: The 
pump photon with energy u)i excites an electron from (a) the 
lower Hubbard band (LHB) [more precisely the Zhang-Rice 
singlet band (ZRB)] or (b) the nonbonding (NB) oxygen band 
to the upper Hubbard band (UHB). The excited electron is 
emitted outside with the kinetic energy i?km by the probe 
photon with energy u>2- 



and the other is due to two sequential (cascade) excita- 
tions, expressed as 
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where N is the number of sites, -Ekin is the kinetic en- 
ergy of photoelectron, and ld\ (102) is the pump- photon 
(probe-photon) energy. The state \i) is the ground state 
of a half- filled system with energy of Ei, and |m) (|/)) 
is the intermediate (final) state with energy of E m (Ef), 
which is the eigenstate of i?Hub- The current operator 
j x connects \i) with |m), where an electric field is as- 
sumed to have a polarization along the x direction in 
the square lattice. We choose the energy relaxation con- 
stant r mm to be twice as large as the phase relaxation 
constant I\ m neglecting the pure dephasing^i The re- 
laxation constants are also assumed to be independent of 
|m): r mm = 2T im = 2T. We chose T = 0.2t. 

In order to calculate I s and I c , we use numerically ex- 
act diagonalization techniques for small clusters with pe- 
riodic boundary conditions. For I s in Eq. @, a standard 
technique combining the conjugate gradient and Lanc- 
zos methods is employed. For I c in Eq. (0, we select 
some dozens of \m) that have large values of (m^li)! 2 . 
The delta functions are broadened by a Lorentzian with 
a width of Q.2t. 

It is important to notice that AR-2PPES spectrum 
generally contains two kinds of energy-dependent spec- 
tral features^ One is the spectrum that satisfies E^ n = 
u>i +u>2+€i, £i being the energy level of an occupied state. 
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FIG. 2: Pump-photon energy u\ dependence of two-photon 
photoemmision spectra at k = (n, 0) in a 4 x 4 t-t'-t"-U clus- 
ter, (a) Simultaneous excitation, and (b) cascade excitation. 
The range of uj\ is determined from the optical conductivity 
cr(u}) in (c). The dashed (dotted) line in (a) and (b) denotes a 
guide to the eyes for the spectrum coming from LHB (UHB). 



Since i?kin is proportional to uj\ + L02 , a peak structure 
obeying this relation is called the 2lo peak, from which we 
can extract knowledge of the occupied state. The other 
gives Skin = <^2 + e m , £m being the energy level of an un- 
occupied state. A peak satisfying this relation is called 
the uj peak and provides knowledge of the unoccupied 
state. 

First, we examine the case where the pump pho- 
ton induces a dipole transition from ZRB to UHB, as 
schematically depicted in Fig. ^a). Since ZRB is re- 
garded as the LHB in the single-band Hubbard model, 
the dipole transition is controlled by a current opera- 
tor obtained from Eq. j x = i ^ k a a^c^ a c^ a with 
ait = —2t sin k x — At' sin k x sin k y + 2t" sin 2k x . 

Figures EJa) and 2(b) show the dependence of I s and 
/ c , respectively, at k = (tt, 0) on the pump-photon energy 
cl>i, calculated by using a half-filled 4x4 cluster where 
fx = 4.55i. We choose the range of uj\ (5t < toi ^ 13t) 
from the optical conductivity <t(uj) shown in Fig. Hfc). 
We note that the range examined corresponds to 1.7 eV ^ 
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FIG. 3: Momentum dependence of two-photon photoemmi- 
sion spectra I a from simultaneous process in a 4 x 4 t-t'-t"-U 
cluster. The pump-photon energy is that (a) w\ = 6t and (b) 
uj\ — 9t. The single-particle spectral function A(k, io) at half 
filling is shown in (c). A(lt,u>) below and above uo/t = can 
be comparable to I s in (a) and (b), respectively. 



to\ ^ 4.6 eV, taking t — 0.35 eV. Since I s and I c are 
plotted as a function of (-Ekin — ^2) A, the 2w peak con- 
taining information on the LHB appears as a structure 
whose position increases linearly with lo\. Such a peak 
is seen along the dashed lines. We find that the 2ui peak 
in I s is clearer than that in I c . This is due to the fact 
that I s includes a virtual excitation via \m) but not in 
I c . As for a structure whose position is independent of 
Wi, we can find a broad peak at around E^ n — L02 = 2.5£ 
in both J s and I c . The weight is enhanced at around 
-Bkin — ^2 = 8.5i, where <j(ui) shows a peak. For the to 
peak, I s and J c show qualitatively similar lu\ dependence. 

Figures E^a) and 3(b) show the detailed momentum 
dependence of the 2io and to peaks in I s at to\ = 6t 
and 9t, respectively. The former, plotted as functions of 
Eki n — uj\— u>2, can be compared to the occupied side of the 
single-particle spectral function j4(k, u>) at half filling 2 in 
Fig. [3f c) , while the latter can be compared to the unoc- 
cupied side. As expected, I s at uj± = 6t shows momentum 
dependence globally similar to that of A(k, to) in the en- 
ergy region of —8 < to/t < 0, although the details of spec- 
tral shape, for example, at k = (0,0) are different. We 



note that the spectral weight below i?kin - = — 8t 

is due to one-photon photoemission keeping an electron 
in the UHB. In I s at i0\ = 9t, despite the broad spectral- 
weight distributions and additional weights at around 
£kin - cj 2 = 3t in the k = (0,0), (tt/2,0), and (0, tt/2) 
spectra, the momentum dependence show a similarity to 
that of A(k, to) for to > 0: From (tt, tt) to (vr, 0) and (0, 7r), 
the spectra show a dispersive feature qualitatively consis- 
tent with the dispersion in the UHB. However, because 
of diffusive features in the spectra, it seems to be difficult 
to identify the location of the bottom of the UHB exactly 
from the co-peak analysis. 

Next we consider the case where the energy of the 
pump photon is tuned to an excitation from the NB 
oxygen band to UHB [see Fig. ^b)]. In order to 
simplify the problem, we construct the Wannier or- 
bitals centered at the copper and oxygen sites ac- 
cording to Ref. 1, and rewrite a current operator be- 
tween copper and oxygen sites by using the Wan- 
nier orbitals. After all, we obtain an expression of 
the cuurent operator between the occupied NB band 
and unoccupied UHB: j x = ^ k a /3k a bk tl7 , where 
6k. a is the annihilation operator of the NB state, and 
/?k = C cos(k x /2) sin(jfej,/2) /y/l — (cosk x + cosk y )/2 
with C = 2d,Q Vl _oTp ( i, dcu-O an d T p d being the distance 
and hopping amplitude between neighboring Cu and O, 
respectively. In the present work, we assume that (i) the 
dispersion of the NB band is negligible^ and (ii) there 
is no interaction between the hole left in the NB band 
and electrons in both LHB and UHB. Under these as- 
sumptions, we can neglect 6k, a m jx- Thus, the dipole 
transition process results in an electron-addition process 
with the momentum-dependent factor /3k- We also take 
the coefficient C to be unity for simplicity. In addition to 
these simplifications, we use the t-t'-t"-J model instead 
of Eq. since only UHB is necessary in this process. 
We note that the t-t'-t"-J model gives almost the same 
results as those of the t-t'-t"-U model when U is large 
enough to satisfy J = At 2 /U. The chemical potential of 
the half-filled t'-t'-t"-J model is determined in order for 
LHB and UHB to be separated by 

In Fig. ^] we show I s obtained by tuning lo\ to an exci- 
tation energy from the localized NB band to the bottom 
of UHB located at (tt, 0) and (0, ir). We note that only the 
(0, 7r) state is occupied by the excitation according to the 
momentum-dependent coefficient /3k- Here lo\ = s p + 2At 
with the level of the NB oxygen band e p . The highest- 
energy structure in / s is located at (0, 7r) and its energy 
is the same as that in A(\c, to). This means that, by tun- 
ing the pump- photon energy to the bottom of the UHB, 
we can observe its position in the momentum and energy 
spaces. If we increase io\ to be tuned to a quasiparticle- 
peak position at k = (47r/5, 2tt/5) with to = 3.7t in 
Fig. EJb), the high-energy edge of the 2P-APRES spec- 
tra appears at (47r/5, 2ir/5) with — L02 — 3.7t (not 
shown). Therefore, the dispersion of the UHB near (tt, 0) 
and (0, 7r) would be detectable by changing o>i, provided 
the dispersion of the NB band is negligible^ We note 
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FIG. 4: (a) Angle-resolved two-photon photoemmision spec- 
tra I a from the simultaneous process in the case that the pump 
photon excite an electron from the NB band, obtained by us- 
ing a \/20 x \/20 t-t'-t"-J cluster. The energy of the pump 
photon is tuned to the energy deference between the NB band 
and the bottom of the UHB with momentum (n, 0) and (0, w) 
as shown in the electron-addition spectral function A(k,Lj) in 
(b). 



that J c shows the same behaviors as those of I a (not 
shown) . 

In addition to the (0,7r) structure, there are spectral 
weights at smaller momentum region with lower kinetic 
energies. For example, a peak appears at (0,0), sepa- 
rated from the (0,7r) peak by ~ t, which is not present 
in A(k, lu). Such additional structures come from the 
following reason. In the process where the NB state is 
excited, the final state |/) belongs to the half- filled sys- 
tem whose low-lying excitations are of the spin wave. 
Actually the (0, 0) peak exists at the eigenstate of the 
spin- wave excitation. Since the peak at (0, tt) comes from 
the ground state of the half-filled system, the energy sep- 
aration between the (0, tt) and (0, 0) peaks is the same 
as the spin- wave width between the momentum transfers 
q = (0,0) and (0,7r). Accordingly the q = (7r,7r) spin- 
wave state contributes to the (tt, 0) peak in Fig. EJa), 
which is thus expected to be degenerate with the (0, tt) 



peak in the thermodynamic limit. After all, we can say 
that the additional states in AR-2PPES contain knowl- 
edge of the spin excitation at half filling. 

Finally we comment on experimental conditions that 
would confirm the present theoretical results. The most 
crucial point is whether one can get -Ekin enough to reach 
to the momentum (tt, 0) and (0,7r). Here we note that 
the maximum of -Ekin is given by E£?^ x = E gap /2 + 
Egap being the Mott-gap magnitude with approximately 
4i according to Fig. Etc). In the case of the excitation 
from LHB to UHB, lo\ should be around 9t from Figs.|3 
andEl If oJ2 = uji, E^ x = lit ~ 4 eV. In the case 
from the NB band, ui\ would be 9i + Ee, where Eb is the 
binding energy of ZRB and approximately 2 eV. Thus, 
by assuming w 2 = Wi, ~ 6 eV in this case. On 

the other hand, the minimum value of Eki n necessary to 
reach (tt, 0) can be estimated to be ~ 6 eV, by taking the 
lattice constant (~ 3 A) and the work function (~ 4 eV) 
into account. This leads to the conclusion that, under 
the condition loi — toi, the momentum (tt, 0) cannot be 
reached for the excitation from LHB but can critically 
for the excitation from the NB band. If we take L02 > <-^i , 
the condition is relaxed and the possibility to observe the 
(tt, 0) state is enhanced. 

In summary, we have proposed that AR-2PPES is a 
promising technique to observe the location of the bot- 
tom of the UHB in 2D insulating cuprates. When the 
pump photon is tuned to an excitation from the LHB 
to UHB, the bottom of the UHB is less clear because of 
diffusive spectral features. On the other hand, when the 
photon energy is tuned to an excitation from NB oxygen 
band, we clearly see the bottom of the UHB. In addition 
to this, additional spectra that are not present in inverse 
photoemission spectra emerge in AR-2PPES. Their ori- 
gin is attributed to the spin excitation expected at half 
filling. 
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